ON THE GLOBAL REGULARITY OF AXISYMMETRIC 
NAVIER-STOKES-BOUSSINESQ SYSTEM 
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Abstract. In this paper we prove a global well-posedness result for tridimen- 
sional Navier-Stokes-Boussinesq system with axisymmetric initial data. This sys- 
tem couples Navier-Stokes equations with a transport equation governing the den- 
sity. 
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1. Introduction 

The purpose of this paper is to study the global well-posedness for three-dimensional 
Boussinesq system in the whole space with axisymmetric initial data. This system 
is described as follows, 



(1) 



dtv + V ■ S/v - Av + Vp = pe^, (t, x) eR+x 
dtp + v - Vp = 0, 
divf = 0, 

V\t=o = Vo, p\t=o = Po- 



Here, the velocity v = {v^,v'^,v^) is a three-component vector field with zero diver- 
gence. The scalar function p denotes the density which is transported by the flow 
and acting for the first equation of ([1]) only in the vertical direction given by e^. The 
pressure p is a scalar function related to the unknowns v and p through an elliptic 
equation. Remark that when the initial density po is identically zero then the above 
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system is reduced to the classical Navier-Stokes system which is widely studied: 

( dfV + V ■ Vf — Av + Vp = 



Recall that the existence of global weak solutions in the energy space for ([2]) goes 
back to J. Leray [23j in the last century. However the uniqueness of these solutions 
is only known in space dimension two. It is also well-known that smooth solutions 
are global in dimension two and for higher dimensions when the data are small in 
some critical spaces; see for instance [21j for more detailed discussions. Although the 
breakdown of smooth solutions with large initial data is still now an open problem 
some partial results are known outside the context of small data. We refer for 
instance to recent papers of J.-Y. Chemin and I. Gallagher [101 IH] where global 
existence in dimension three is established for special structure of initial data which 
are not small in any critical space. 

There is an interesting case of global existence for ([2]) corresponding to large initial 
data but with special geometry, called axisymmetric without swirl. Before going 
further in the details let us give some general statements about Navier-Stokes system 
in space dimension three. First, we start with introducing the vorticity which is a 
physical quantity that plays a significant role in the theory of global existence; for a 
given vector field v the vorticity u is the vector defined hj u = curlf. Thus we get 
from the vorticity equation 



According to Beale-Kato-Majda criterion the formation of singularities in finite time 
is due to the accumulation of the vorticity. In other words, to have global existence it 
suffices to bound for every time the quantity || However the main difficulty 

arising in dimension three is the lack of information about the manner that the 
vortex-stretching term u ■ V v affects the dynamic of the fluid. 
For the geometry of axisymmetric flows without swirl we have a cancellation in the 
stretching term giving rise to new conservation laws. We say that a vector field v is 
axisymmetric if it has the form: 

v{t,x) = v^{t,r, z)er + v^{t,r, z)ez, x={xi,X2,z), r = {xf + xl)^ ^ 

where (e,-, eg, e^) is the cylindrical basis of and the components and do not 
depend on the angular variable. The main feature of axisymmetric flows arises in 
the vorticity which takes the form. 



(2) 




dtuj + V ■ Vu; — Au; = u ■ Vw. 



UJ = {dzV^ — drV^)eg 



and satisfles 



dtUj + V ■ Vuj — Auj 
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Since the Laplacian operator has the form A = drr + ^dr + dzz in the cyhndrical 
coordinates then the component ujq of the vorticity will satisfy 

(3) dtue + V ■ Vuje - Acu^ + ^ = -^^e- 

Consequently, the quantity T := ^ obeys to the equation 

dtT + vVT-AT- -drT = 0. 

r 

Obviously, we have for p > I, J ^drT\T\P~^ {sign T)dx < and then we deduce that 
for all p E [l,oo] 

l|r(t)|Up < iifoIUp. 

It was shown by M. Ukhoviskii and V. Yudovich [25j and independently by O. A. 
Ladyzhenskaya [20] that these conservation laws are strong enough to prevent the 
formation of singularities in finite time for axisymmetric flows. More precisely, the 
system ([2]) has a unique global solution for vq G such that ujq,^ G Pi L°°. 
Remark that in term of Sobolev regularities these assumptions are satisfied when 
Vq G H'^ with s > |. Few decades afterwards, S. Leonardi, J. Malek, J. Necas and M. 
Pokorny [22] weakened the initial regularity for vq G H^. This result was recently 
improved by H. Abidi [1] for vq G H2. 

Let us now come back to our first problematic which is the study of global well- 
posedness for the Boussinesq system ([T]). In space dimension two many papers are 
recently devoted to this problem and the study seems to be in a satisfactory state. 
More precisely we have global existence in different function spaces and for different 
viscosities, we refer for example to[2l[7l[Hl[13l[Tl[T5l[l6l[I7l[T8l[l9]. 
In the case of space dimension three few results are known about global existence. 
We recall the result of R. Danchin and M. Paicu [H]. They proved a global well- 
posedness result for small initial data belonging to some critical Lorentz spaces. 
Our goal here is to study the global existence for the system ([1]) with axisymmetric 
initial data, which means that the velocity vq is assumed to be an axisymmetric 
vector field without swirl and the density po depends only on {r,z). It should be 
mentioned that this structure is preserved for strong solutions in their lifespan. 
Before stating our main result we denote by II^ the orthogonal projector over the 
axis [Oz). Our result reads as follows. 

Theorem 1.1. Let vq G be an axisymmetric vector field with zero divergence 
and such that ^ G L^. Let po E L'^ n L°° depending only on {r, z) and such that 
supp Po does not intersect the axis (Oz) and Ilz{supp po) is a compact set. Then the 
system ([1]) has a unique global solution {v,p) such that 

veC{R+-H')nLl^ (M+;iy^--), 

-eL^M+:L'), peL^M+;L'nL^). 
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The assumption that the density p is zero in some region of the space M'^ is not very 
meaningful from a physical point of view. However we can relax this hypothesis and 
extend our result to more general case: the initial density is assumed to be constant 
near the axis {Oz) and for large value of z. More precisely, we have the following 
result. 

Corollary 1.2. Let vq G he an axisymmetric vector field with zero divergence 
and such that ^ G L^. Let po ^ L°° depending only on (r, z) such that po = Cq, 
for some constant cq in a region of type |x; r < tq, \z\ > |-2o|}, with tq > and 
Po ~ Co & L"^ . Then the system ([T]) has a unique global solution (f,p) such that 

vECiR+;H')nLl^ iR+;W'n, 

- G (M+; L^), p - Co G (M+; L' n L^). 
r 

The proof of this corollary is an immediate consequence of Theorem 1 Indeed, we 
set p(t,x) = p(t,x) — Co, then the system ([1]) is reduced to 

{dtv + V ■ Vf — Av + Vp — coCz = pcz 
dtp + v-Vp = 
divv = 
v\t=o = vo, p\t=o = Po- Co- 

Now by changing the pressure p to p = p — coz v/e get the same system ([1]) and 
therefore we can apply the results of Theorem 11.11 

Next we shall briefly discuss the new difficulties that one should deal with compared 
to the system ([2]). First we start with writing the analogous equation to for the 
vorticity. An easy computation gives 



curl(pe^) = 




= -{drp)e0. 



This yields to 

(4) dtujg + V ■ Vug - AuJe + = —^e - drp. 

It follows that the evolution of the quantity F := ^ is governed by the equation 

(5) dtT + vVT-AT-'^drT = -^- 

Since the density p satisfies a transport equation then the only conserved quantities 
that one should use are ||p(t)||LP for every p G [l,C)o]. Loosely speaking, the source 
term ^ in the equation that one need to estimate has the scale of Ap and 
then to estimate F in some Lebesgue spaces one can try for example the maximal 
regularity of the heat semigroup. However it is not at all clear whether we can prove 
a suitable maximal regularity because the involved elliptic operator —A — -dr has 
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singular coefficients . Now, by taking the L^-inner product of with F we are led 
to estimate the quantity ||p/r||j;^2 which has the scaling of ||Vp||i2. In other words 
we have 

l|r(t)|li.+ ri|r(r)ii5,,rfr< ||roiii.+ f \\ipM{r)\\hdr. 

Jo Jo 
Remark that the quantity ||p/r||j;^2 is not well-defined if we don't assume at least 
that p{t, 0, z) = 0. As we have seen previously in the statement of Theorem 11.11 we 
need more than this latter condition: the initial density po has to be supported far 
from the axis (Oz). To follow the evolution of the quantity ||p(i(:)/r||/^2 the approach 
that consists in writing the equation of ^ fails because it gives an exponential 
growth, 

||p(t)/r|U^<||po/r|U2e"^'/''"-J--. 

Our idea to perform this growth relies in studying some dynamic aspects of the 
support of p{t) which is nothing but the transported of the initial support by the 
flow. We need particularly to get a proper low bound about the distance from 
the axis (Oz) of the support of p(t), see Proposition 13. 2[ To reach this target and 
once again the axisymmetry property of the flow will play a crucial role since the 
trajectories of the particles are contained in the meridional plane. This approach 
allows us to improve the exponential growth for ||p(t)/r||j;^2 to a quadratic estimate. 
Roughly speaking we obtain the estimate 

\\pit)/r\\L^ < Co\\v''/r\\LiL-°{^ + \Ml}l--)- 

The rest of the paper is organized as follows. In section 2 we recall some ba- 
sic ingredients of Littlwood-Paley theory. In section 3 we study some qualitative 
and analytic properties of the flow associated to an axisymmetric vector fleld. In 
section 4 we give some global a priori estimates. The proof of Theorem 11.11 is done 
in Section 5. 

2. Preliminaries 

Throughout this paper, C stands for some real positive constant which may be 
different in each occurrence and Co for a positive constant depending on the initial 
data. We shall sometimes alternatively use the notation X <Y for an inequality of 
type X < CY. 

Let us start with a classical dyadic decomposition of the whole space (see [9]): there 
exist two radial functions x ^ 'P(K^) and ip G I?(M'^\{0}) such that 

(1) x(0 + E^(2-'^^) = l V^GM^, 

q>o 

(2) supp </?(2-P-) n supp ^(2-'?-) = 0, if \p-q\> 2, 

(3) g > 1 ^ suppx n supp v^(2^^) = 0. 
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For every u G iS'(]R^) we define the nonhomogeneous Littlewood-Paley operators by, 
A_iM = x(P)u; Vg G N, AqU = f{2~m)u and SgU = ^ Aju. 

-l<j<q-l 

One can easily prove that for every tempered distribution u, 



qU. 



(6) « = E ^ 

In the sequel we will make an extensive use of Bernstein inequalities (see for 
example [9j). 

Lemma 2.1. There exists a constant C such that for kEN,l<a<b and u G 
we have 

sup iia^viU" < C^2'?('=+=^(^-^))||VlU''. 

\a\=k 

and for g G N 

C~'=25^'||Agu||La < sup Wd'^AquWLa < CH'^^AquWic^. 

\a\=k 

Let us now introduce the basic tool of the paradifferential calculus which is Bony's 
decomposition It distinguishes in a product uv three parts as follows: 

uv = TuV + Tr^u + IZiu, v), 

where 



TuV = Sq^iuAgV, and 7l{u, v) = AquA 



qV, 



1 



with Ag = E 

i=-l 

Ty_v is called paraproduct of f by m and Tl{u, v) the remainder term. 

Let (p, r) G [1, +cxd]^ and s G M, then the nonhomogeneous Besov space 5*,. is the 

set of tempered distributions u such that 



\u\\bs^ := (2'^'\\Aqu\\LP^ < +00. 



We remark that the usual Sobolev space H'^ agrees with Besov space -B2 2- Also, by 
using the Bernstein inequalities we get easily 

s+3(-^-^) 

Bp^^^^ "-^ Bp2,r2''^ , Pi<P2 and ri < r2. 
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3. Study of the flow map 

The main goal of this section is to study some geometric and analytic properties of 
the generalized flow map associated to an axisymmetric vector field. 

ip(t,s,x)=x + J v{t,iP{t, s,x))dT. 

This part will be the cornerstone of the proof of Theorem II. 1[ We need first to 
recall some basic results about the generalized flow. If the vector field v belongs to 
Ljjjj,(]R, Cl) then the generalized flow is uniquely determined and exists globally in 
time. Here denotes the space of functions with continuous bounded gradient. 
In addition, for every t, s G M, ip{t, s) is a diffeomorphism that preserves Lebesgue 
measure when div v = and 

ip~'^{t, s, x) = ip{s, t, x). 
Now we define the distance from a given point a; to a subset A C by 

d(x, A) := inf ||x — y\\, 

where || ■ || is the usual Euclidian norm. The distance between two subsets A and B 
of is defined by 

d(A,B):= inf lla; — -ull. 

The diameter of a bounded subset A C is defined by 

diam A = sup — ?/||. 

Our first result is the following. 

Proposition 3.1. Let v be a smooth axisymmetric vector field and ip(t,s) its flow. 
Let X ^ {Oz) and r{x) := d{x, (Oz)), then 

(1) For every s G M, the trajectory r^,. ,, := \^ilj{t, s,x),t G M} is a smooth curve 
contained in the meridional plan. 

(2) For every s G M, the trajectory T^^s does not intersect the axis (Oz), that is 
^x,s n (Oz) = 0. More precisely, 

r{x)e-\^s\\'T(^)U^dr\ < d{ip{t,s,x),{0z)) < r(x)el^''*ll'f(")l'^°°'^"'. 

Proof. (1) We start with the decomposition of the vector x in the cylindrical basis: 

^cos Ox'^ 

X = r{x) I sin6'2; | + ZxCz, with rx > 0. 
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We decompose also the generalized flow ilj{t, s,x) in the cylindrical coordinates 

(cos{9{t, s, x))\ 
sm{6{t,s,x)) + z{t,s,x)ez. 

Since v is smooth in space-time variables then the flow map ip{t, s, x) is also smooth. 
Now we intend to prove that r(t,s,x) remains strictly positive for all t,s e R. 
Assume that there exist ti, Si G M such that ip{ti, Si, x) belongs to the axis {Oz). As 
the restriction of the vector fleld v on the axis (Oz) satisfles v{0, 0, z) — v^{t, 0, z)ez 
then trajectory for every xq G {Oz) lies in this same axis 



ip{t, s, xo) = xo + (O, 0, v^(r, s, Xo)dT) . 



Clearly one can choose Xq such that iplti, Si,Xq) = ilj{ti,Si,x) and this contradicts 
the fact that the flow map is an homeomorphism. Consequently and from the 
smoothness of the generalized flow map one can prove easily that the functions 
r{t, s, x), 9{t, s, x) and z{t, s, x) are also smooth in each variable. This allows us to 
justify the following computation 

(cos{9{t, s, x)y 
sm(9{t, s, x)) 

/-sm{9{t,s,x))\ 
+ r{t, s,x)dt9{t, s,x) I cos{9{t, s,x)) I + dtz{t, s,x)ez. 

\ / 

Since the vector fleld v is axisymmetric without swirl then 

/cos(^(t,s))\ 

s, x)) = r(t, s, a;), s, x)) sin(6'(t, s)) + v^{t,r{t,s,x),z{t,s,x))ez. 

\ / 

Thus we get by identiflcation 

dtr{t, s, x) = v''{t, r{t, s, x), z{t, s, x)), 
(7) r{t,s,x)dt9{t,s) =0, 

dtz{t, s, x) — v^{t, r{t, s, x), z{t, s, x)). 

Prom the above discussion we have r{t, s,x) > for all s G IR and therefore we get 

9{t, s, x) = 9{s, s, x) = 9^, Vi, s G R. 

It follows that for every s G R the trajectory T^^s '■= {'V'(^) •s, a;), i G R} lies in the 
meridional plan. 
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(2) From the first equation of (171) we get 

r{t, s, x) = r{x) + / r{T, s, x)~^v^ (^T,r{T, s), z{t, s)')r{T, s, x)dT. 

J s 

Using Gronwall lemma we get 

r{t,s,x) < r(x)el^»*ll'f 

Tliis gives tlie second inequality of the r.h.s since di^ipit, s,x), (Oz)) = r(t,s,x). 
Now we apply the above inequality by taking ip{s,t,x) instead of x. Since 

ip(t, s, ip{s, t, x)) = X then 

Interchanging (s,t) and {t,s) gives 

^(a;)e-l/.*ll^«lli-'^^l < r{t,s,x). 

This achieves the proof of Proposition 13.11 □ 

Proposition 13.11 will be of much use in deriving some a priori estimates about solu- 
tions of transport equations. The first application is given below, 

Proposition 3.2. Let v be a smooth axisymmetric vector field and p a solution of 
the transport equation 

( dtp + v-Vp = 
\ P\t=o = Po- 

(1) Assume that (i(supp po, (Oz)) = ro > 0. Then we have for every t > 

d(supp (Oz)) > ro e-^o\\i-(r)\\Loodr_ 

(2) Denote by Hz the orthogonal projector over the axis (Oz). We assume 
that n2(supp Po) is compact set with diameter d^. Then for every t > 0, 
n2(supp p{t)) is compact set with diameter d{t) such that 



d{t) <do + 2 [ 
Jo 



Proof. (1) The solution p of the transport equation is completely described through 
the usual flow ip, that is, p{t,x) = po{'ip~^{t,x)). Here ip{t,x) := ip{t,0,x) where 
i(j{t,s,x) is the generalized flow introduced in the beginning of this section. Thus 
it follows that supp p{t) = ipit, supp po). Let now y G supp p(t) then by deflnition 
y = ipit, x) with X e supp po- Therefore d{y, (Oz)) = rit, x) with r(t, x) := r{t, 0, x) 
and so from Proposition 13.11 we get 

d{y,{Oz)) > d{x,{Oz))e-^o\\'^i-)\\L^<ir 

> C/(SUPP Po, (Oz)) e-^O ll^Wboodr 

> roe-^oll^MIU--^-. 
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This concludes the first part. 

(2) Let x,x & supp po and denote by y{t) = ip(t,x) and y{t) = ip{t,x) . We set 
successively z(t) and z(t) the last component of y(t) and y(t). Using the equation 
(I7j) with s = we get 

i(t) = v'{t,r{t,x),z{t)). 
Integrating this differential equation we get 

^(t) = ^(0) + f v''{T,r{T,x),z{T))dT. 
Jo 

It follows that 

\z{t) - Z{t)\ < \Z{0) - 5(0)1 +2 [ \\v{T)\\LoodT. 



This yields 

diam(n2(supp p(t))) < diam(n2(supp po)) + 2 / || dr. 

Jo 

The proof is now completed. □ 

Now we are in a position to give an estimate of the quantity ||p(t)/r||/,2 which is very 
crucial in the proof of Theorem 11.11 The growth that we will establish is quadratic 
and this improves the exponential growth that one can easily obtain by writing the 
equation of p/r. More precisely we have 

Corollary 3.3. Let v be a smooth axisymmetric vector field with zero divergence, 
Po E L?' r\ L°° and p be a solution of the transport equation 

( dtp + v-Vp = 
\ p\t=o = Po- 

Assume in addition that 

(i(supp Po, (Oz)) := To > and diam(n2(supp po)) := rfo < oo- 

Then we have 
p\t,x) 



Ux<\\\pYl^+MPo\\U I \\{v'-/r){r)\\L^dT(do + 2 f \\v{T\\Loodr), 
fo Jo ^ Jo ^ 



with r = [x\ + x\) 2 . 
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Proof. We have from the definition 



J r>ro ' J r<rx) ' 

< illpWIIi^ + llpWIIio. / ^dx 

'0 J {r <ro}nsupp p(t) ' 

1 „ „2 „ u2 f 1 , 

< — IIPo||l2 + IIpoIIl- / —dx. 

We have used the conservation of the L°°-norm of p. Now using Proposition 13.21 we 

get 

/ \dx < 27r( / ^ -dr)( I dz 

< 27r^ ||y(r)||ioo(ir(^do + 2^ ||t^(T||Loodr' 
This concludes the proof. □ 



4. A PRIORI ESTIMATES 

This section is devoted to the a priori estimates needed for the proof of Theorem 1 
We distinguish especially two kinds: the first one deals with some easy estimates 
that one can obtained by energy estimates. However the second one is concerned 
with some strong estimates which are the heart of the proof of our main result. 

4.1. Weak a priori estimates. We will prove the following energy estimates. 

Proposition 4.1. Let vq & L'^ be a vector field with zero divergence and po G 
n L°°. Then every smooth solution of ([T]) satisfies 

\\p{t)\\L^nL^ < ||Po||L2nL°°- 

\\vmh+ f\\Vv{T)\\Ur<Co{l+t'). 
Jo 

Proof. The first estimate is obvious since the fiow preserves Lebesgue measure. For 
the second one we take the L^-inner product of the velocity equation with v. Then 
we get after some integration by parts 

(8) ^|ll^^Wlli. + l|V^(t)||L<l|t;(t)lU.lb(t)IU.. 

After simplification that one can rigorously justify, the last inequality leads to 
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Integrating in time this inequality yields 

HmL^ < ho\\L^+ [ \\p{r)\\L^dr. 
Jo 

Since ||p(t)||L2 = ||po||l2, then 

||t'(t)||L2 < ||^^o||l2 + t||po||L2- 

Putting this estimate into ([8]) gives 

^II^WIli2 + ^ \\Vv{T)\\l,dT < ^WvoWl. + {\\vo\\L^+t\\po\\L^)\\po\\L^t. 

This gives the desired estimate and the demonstration of the proposition is now 
accomplished. □ 

The next proposition describes some estimates linking the velocity to the vorticity 
by the use of the so-called Bio-Savart law. 

Proposition 4.2. Let v be a smooth axisymmetric vector field with zero divergence 
and denote ui = ujgCg its curl. Then 

(1) 

1 1 
\\v\\lo^ < C\\iUe\\'l2\\uJe\\^^i. 

(2) 

1 i_ 

\\v''/r\\L^ < C\\uJe/r\\l2\\ujg/r\\^^^. 

Proof. To start with, recall the classical Biot-Savart law 
It follows that 

Hx)\ < 

Let A > be a real number that will be fixed later. Decompose the convolution 
integral into two parts as follows 



x-y\<x F y\ A^-y\>^ 1^ y\ 

Jl{x) + J2{x). 
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Now using Holder inequalities we get 



|Ji||loo < llcuellisf / -^(ix 



\x\<X \x\ 



5 



For the second integral we use again Holder inequalities 



|J2|U°° < lke||L2( / T-udx 



1 ■ ' 



>A Fl 



This yields 

By choosing A = l^"!!^^ we get 

1 1 

It suffices now to use Sobolev embedding H^(R^) ^ to get the desired result. 
For the second estimate we recall the following inequality, see for instance [51 [21] , 

Iv'^/rlix) < -k \uJe/r\){x). 

Now we conclude similarly to the first estimate. □ 

4.2. Strong a priori estimates. The task is now to find some global estimates 
about more strong regularities of the solutions of (Q. The estimates developed 
below will be the basic ingredient of the proof of Theorem 11.11 

Proposition 4.3. Let Vq G be an axisymmetric vector field with zero divergence 
such that ^ e L'^ D L°° depending only on {r, z) such that supp po 

does not intersect the axis (Oz) and Ilz{supp po) is a compact set. Then every 
smooth solution {v, p) of the system ([T]) satisfies for every t > 0, 



JO 

Moreover, for every p g]3, 00] we have 

(9) \\v\\ + WVvhr^^ < Coe^^P^«*' 

r 1 nP t 

The constant Co depends on the initial data. 
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Proof. Taking the L^-inner product of the equation (jll) with ug we get 

^ -rrll'^e II ^2 + II Vt^e 11^2 + 11 — 11^2 = / v''—ujgdx~ [ drpuedx. 
2 at r J r J 

For the first integral term of the r.h.s we use Holder inequalities 

r^e J , 

V — ujedx < ||'y||Le|| — ||l2 ||u;e||2,3. 

For the second term we integrate by parts taking into account that the density p 
vanishes on the axis (Oz) accordingly to Corollary (13.31) 

— J drpojgdx = —271 J OrPi^e^dr dz 

= 2t[ j pdrtoerdr dz + 271 J p—rdrdz 

= J p[drUJe + '^)dx. 

It follows from Holder inequalities 

- j drpUJgdx< ||p||i2 (11^11^2 + ||(9^C^0||i2). 

Putting together these estimates yields 

^ II l|2 , IIV7 i|2 , ||^''||2 ^ II II ll^eil II II 

2^ll^e||L2 + IIV^elL^ + ||ylL2 < ll^llL6||yL2ll^e||L3 

fw^eii , 11^ II 

+ ll/'ll^^l^ll^llLa + l|c'rt^e||L2 

Using the inequality |a6| < + ^b"^ for the last term we obtain 

rW^^W , lia II \ ^ II l|2 I , l||rj ||2 

nllylL^ + ll^-^elUO < \\p\\L^ + ^\\ — \\L2 + ^\\drUJ9\\L2 



IIpIIl^i 



< llplli^ 



i^ + ^llyli:2 + i||vc..|ii2. 



We have used in the last line the fact that in the cylindrical coordinates 

||Vti;6»||i2 = ||(9rt^6»||i2 + ||(9za;6»||i2. 

Therefore we get 

(10) j^Mh + l|V^.||i2 + ll^ll', < 2\\vhe\\'^\\J\u;ey. + 2\\p\\l,. 



ON THE AXISYMMETRIC BOUSSINESQ SYSTEM 15 

1 i_ 

Combining the interpolation estimate Ht^elU^ ^ Ht^ell^a II Vc<j6i||22 with Young inequal- 
ity \ab\ < Crja^ + ife^, Wi] e]0, 1[ leads to 

ML^WylUl^oWv^ < C|l^llillylliHlli + ^l|Vc.,||i. 
(11) < C'||^,||i. + i||V..,||i. + C||t;||i.||i^||^,. 

Inserting (fTTi) into (fTUI) and using the estimate ||p(t) 11^,2 < UpoIIl^ 

|||^,||i. + i||V^,||i. + ll^ll^, < Cll^.lli. + C||.;||ie||^||'. + 2||po||i2. 
Using Gronwall inequality and Proposition 14.11 we get 

\\Mmh + \\^^o\\%L^ + \\^\\l.,, < Ce^*(||po|li.+ ril^(r)|lia||^(r)||^,cir 
t r t ^ Jo ^ 

< C^e^*^IIPo|li. + ||^^|lW.||^||U. 



r 



L2 



(12) < Coe^'(l + \\-\\l 

Let us now move to the equation and try to estimate in a proper way the 
quantity F := Then taking the L^-inner product of ([5]) with F and integrating 
by parts using the incompressibility of the flow and p{t, 0,z) = gives 

^^llrlli. + ||9,F||i. + ||9,F||i. -An J dr{T)Tdrdz = -2n J drpTdrdz 

2n I -drVrdrdz 
r 



< ||-|L.I|5.r|U.. 



Since 



//• r+oo 
dr{T)Tdrdz = 2tx j j dr{Tfdrdz < 0, 
Jr Jo 

then and once again by Young inequality yields 

^ll^^^ll^ -I- II Vf— ) 11^ < 11-11^ 

^^W ^ II V / IIl2 — II ^ IIl2 

Integrating in time this differential inequality we get 

(13) II^WIll. + /' l|V(^) W||> < ||^(0)||:, + /' \\l(r)\\lir. 

r Jo ^ " Jo ^ 

To estimate the last term of the above inequality we will use Corollary 13.31 



|^(t)||^,<Co + Co^*||^(r)||,^rfr(l + ^* 



\v{T)\\LoodT]. 
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Thus 

ll?W|lL + l|v(^)||l,.,<c„(i + «) 

-t , l-t' 



+ Co^{^ ||(^)(r)||Lo^rfr(l + ^ ||i;(r)|Uocrfr)}cit' 

<Co(l + t) + Cot r||(-)(r)|Uoorfr 

+ Cot [ \\i-)ir)\\Lo^dr [ \\viT)\\L^dT. 
Jo ^ Jo 

From Young and Holder inequalities we get successively 

Cot f U-){T)h^dT < [\cot'+\\{-){r)\\l^)dT 
Jo ' Jo " 

< Cot'+ f \\{-){r)\\l^dT 
Jo r 

and 

Cot\\v''/r\\LlL^\MLlL-- < Cot^\\v'/r\\L2L^\\v\\L2Loo 



< Cot' f\\C-^){T)\\l^dT+ f\HT)\\l^dT. 







Putting together these estimates yields 

(14) + /'||t'(r)l|i„<iT, 

Jo 

Recall from Proposition 14.21 that 

W" / ' ~ II ^ IIl2|| V I, ^ )\\^2- 
Accordingly, we get by Young inequalities 

Co{l+t') f\\{v^/r){r)\\l^dr<Co{l+t') f \\-{r)\\ldr+l T ||V(-)(^ 
Jo Jo r 2 Jo r 

Inserting this estimate into ( |T4l) 

II^Wr.. + l|vOl|^.,<C„(l + *»)/||^(r)||> 

(15) + / \\v{T)\\l„dT. 
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Now using Gronwall inequality we find 

(16) + W'^OWln- s Co'*'(i + /' IK'(-)lli-'^-). 

PuttingdU]) into yields 

(17) WMrnh + llV^.lli.^. + < Coe^"*' (l + HT)\\l^dr). 

To estimate the last term in the r.h.s of (ITTl) we use Proposition 14.21 

1 1 
h\\L°- < ll^e|l22||Va;6i||22- 

It follows 

(18) f \\v{T\\l^dT< [ \\uJe{T\\L4^UJg{T)\\L2dT. 

Jo Jo 
Therefore we get from ( fT71) and Young inequality 

\\ooeml2 + \\VuJe\\l2L2 + \\ — \\l,^, < Coe^"*Yl + / ||cue(r|U2 1| Vcu,(r) lU^t^r) 
t r t ^ Jo ^ 

It suffices now to use Gronwall inequality 



1 1 2 

Since in cylindrical coordinates we have ||Vco'||^2 = ||Vco'6i||^2 + Hv'ILa then 

(19) ||o.(^)||i2 + ||Va;||i2,2<Coe-P^o*^ 
Inserting this estimate into (ITS!) gives 

(20) ||t;|L2^o. <Coe^^P^°*'. 
Combining this estimate with f|T6l) and f|T9l) yields 

II^W|i:. + ||v(^)W||i,,.<c„.-«'. 

This concludes the first part of Proposition 14. 3[ Let us now show how to prove the 
estimate ([9]). Let g G N and set Vg := AgV. Then localizing in frequency the first 
equation of ([1]) and using Duhamel formula we get 



(t) =e%(0)+ [\^'-^^''Ag{V{vVv)iT)dT+ re(*--)^A,(P(pe,))(r)t/r, 
Jo Jo 
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where V denotes Leray's projector over solenoidal vector fields. Now we will use a 
local version of the smoothing effects of the heat semigroup, for the proof see for 
example [9], 

||e*^A,P/|Up < Ce-'^*2'lAJ|U., Vp G [l,oo]. 
Combining this estimate with Bernstein inequality gives 



e-'^'-^>''\\A,piT)\\L.dT. 







Integrating in time and using convolution inequalities lead to 

'I^JIlJlp <2-''^ll^a(0)||LP + 2-« / \\A^{v®v){T)\\Lr>dT + 2-^'' [ || A,p(r) |Uprfr. 

^0 Jo 

It follows that 



l^^ll ^ 1+1 < ||A_if ll^iip + ||t;o|| i_i + [ \\{v^v){t)\\ idr 
^Ki <i Jo <i 

*||p(r)|| |-,rfr. 

<i 

For the first term of the r.h.s we combine Bernstein inequality with Proposition 14.11 
(P>2) 

For the last term of the r.h.s we use the embedding ^ B^^^ , for p > 3, combined 
with Proposition 14.11 

t 

\\p{r)\\ i^_,dT < t\\p\\L^Lp 

Bp,i 

< Cot. 

On the other hand we have from Besov embedding 

II^OII 3 1 II'^OII 1 

^p,l "2,1 

< W^oWm- 
Putting together these inequalities we find 

\\vl ,^,<Coil + t')+ I \\{v v){t)\\ dr. 
^t^p,i Jo 



LIB!;, ■ ■ ./n 
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Now we use Besov embeddings, law products and interpolation results 

||t> (g) t>|| 3 < 111; ® f II 3 



< 

< 



f room 3 
I \\L II 11^^, 



2,1 



B? 



2.1 



^ II^IU°=||^^||l2 + ||f ||L-'||t^||f2||a;P 



!_ 1 
2 II, ,11 2 

lL2|l'^ll//i- 

This yields according to Holder inequalities 

1 



1 1 

2 



P,i 



1 1 

2 



It suffices now to use ( IT9l) and (120 



L^°°L2||',^||j;^2jyl- 



Hence we obtain 

\\v\\ 3+1 < Coe^^P^o*'. 

The last estimate of Proposition ([9]) is a direct consequence of the 
embedding B^^ W^''^ and the proof is then accomplished. □ 



5. Proof of the main result 

The existence part can be done in a classical way for example by smoothing out the 
initial data. However we have to use an appropriate approximation that does not 
alter never the initial geometric structure nor the uniform estimates in the space of 
initial data. We will work with the approximation of the identity and show that it 
has the requested properties: let be a smooth positive radial function with support 
contained in 5(0, 1) and such that = 1 in a neighborhood of zero. We assume 

that / (f)(x)dx = 1. For every n G N* we set 0ri(x) = n^(j){nx) and we define the 

family 

Vo,n = 4>n-*^Vo and po,n = 0n * PO- 

We will start with the following stability results. 

Lemma 5.1. (1) Let Vq G be an axisymmetric vector field with zero diver- 
gence and such that (curl Vo)/r G L^. Then for every n eN* the vector field 
fo,n ^-5 axisymmetric with zero divergence. Moreover, there exists a constant 
C depending on cf) such that 

hoAn^ < ll^ol^i, ||(curl VQ^n)/r\\^:, < C||(curl VQ)/r\\^^ 
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(2) Let V be a smooth axisymmetric vector field with zero divergence and p be a 
solution of the transport equation 

dtp"" + V ■ Vp" = 



P \t=0 — PO,n- 

Assume that d{supp po, (Oz)) := ro > and diam (nj;(supp po)) := do < oo. 
Then there exists constants hq and C depending only on r^^do and such 
that 

f (p")^(t,x) 1 

/ — dx < — ||po||i2 +2||po||L-(ln2+ IK'V'^IIlil-) (c^o + II^^IIlil-)- 

Jr3 r Tq 

Remark 5.2. The estimate of the part (2) of the above lemma is httle bit different 
from part (2) of Corollary 13.31 because we have an additional linear term ||f H^^i^^oo. 
Nevertheless, the presence of this term does not deeply affect the calculus seen before 
in the a priori estimates; we obtain at the end the same estimates of Proposition 



Proof. (1) The fact that the vector field fo,„ is axisymmetric is due to the radial 
property of the functions 0„, for more details see [3]. The estimate of vo,n in is 
easy to obtain by using the classical properties of the convolution operation combined 
with ||0n||Li = 1. Concerning the demonstration of the second estimate, it is more 
subtle and we refer to |5] where it is proven for more general framework Lebesgue 
space, that is L^, for all p E [l,oo]. 

(2) In order to establish the desired estimate we need to check that the estimates 
of Proposition 13.21 are stable with respect to n. More precisely, we will prove that 
for sufficiently large n > uq 

(21) ci(supp p"(t),(Oz)) > !°e-^oll^Mlli-'^- 
and 

(22) rf„(t) < 2^0 + 2 /" \\v{T\\LoodT with rf„(t) := diam (n^(supp p"(t))). 



Assume that we have the inequalities (1211) and (1221) then reproducing the proof of 
Corollary 13.31 we get 

1 f 1 

\\ipn/r){t)\\l2 < — llPo.nllia + ||po,n||i- / —dx 
1 „ 1,2 „ 1,2 /" 1 , 

< — IIPo|Il2 + IIPoIIloo / —dx 

'"O J{r<ro}nsupp p„(t) ^ 



< — 

r, 



1 /"* 

-2llPo|li2 + ||po|lioo(ln(2) + ||t;7r||^i^oo)(2c/o + 2 / ||t;(r||Looc/r) . 

"o Jo 
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This is what we want . Now let us come back to the proof of the inequahty ( 1211) . Fol- 
lowing the same proof of Proposition 13.21 it appears that one needs only to establish 
the stability conditions of the support, that is, for large n > uq, 

(23) d(suppp^,(Oz))> ^. 

From the definition we have 

supp Po C ja;; Us < ^ for some , \\y'\\ > roj, with y = {y\y^),y' G M^. 

It is obvious that for x G supp Pq we have — < ^ and then by the triangular 
inequality 

„ /„ 1 

\\x > ro 

Thus there exists ng G N such that for every n > no we have > y- This proves 
that supp Po C {x, ||a;'|| > y} and then we deduce f l23l) . 

Let us now show how to get the estimate ( l22l) . According to part (2) of 
Proposition 13.21 we have 

dn{t) < do,n + 2 / \\v{T)\\Lo^dT, 

Jo 

with (io,n := diam (n2(supp Po,n))- Now it is easy to see from the support properties 
of the convolution that 

2 

do,n < do -\ 

n 

Taking n large enough, n> rii for some ni, yields (io,n < 2(io- Now to be sure that 
both inequalities (|2T1) and (|22l) are simultaneously satisfied we take n > max{no, rii}. 

□ 



Let us now come back to the proof of the existence part of Theorem ll.il We have just 
seen in Lemma 15.1! that the initial structure of axisymmetry is preserved for every 
n and the involved norms are uniformly controlled with respect to this parameter 
n. Thus we can construct locally in time a unique solution 6n) that does not 
blow up in finite time since the Lipschitz norm of the velocity is well controlled as it 
was stated in Proposition 14.31 By standard arguments we can show that this family 
converges to {v, 6) which satisfies in turn our IVP. We omit here the details which 
are very classical and we will next focus on the uniqueness part. Set 
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Let (f*,p*) G Xt, 1 < i < 2 be two solutions of the system ([T]) with the same initial 
data (fo, 9o) and denote 6v = v'^ — , 69 = 9"^ — 9^ . Then 

{dfSv + f ^ ■ V6v — A6v + V6p = —6v ■ Vv^ + Spe^ 
dt6p + v'^ ■ V6p = -6v ■ Vpi 
divt>* = 0. 

Taking the L^-inner product of the first equation with Sv and integrating by parts 
- — 11^2 + II V(5t'|||2 = — y 5v-'Vv^Svdx + J Spe^Svcix 

< \\6v\\l4v'\\l^\\V6v\\l2 + \\6p\\H-AM\m 

< {\\Svh4v'h^ + \\6p\\H-^)\\vsvh. + \\6py-4H\L^- 

Applying Young inequality leads to 

(25) j^\\6v\\l, + ||V5^;||i2 < ||5^;||i.(||t;1i. + 1) + \\Sp\\l-.. 

To estimate ||5p||//-i we will use Proposition 3.1 of [4J: for every p E [2, oo[ 
||5p(t)lk-i < C\\Sv ■ Vp'hiH-^ exp iC\\Vi 

As div 6v = 0, then 

< IIpoI|l-'||5^||lil2. 

Inserting this estimate into fl25l) . we obtain 

|llMi)lli. + l|VMi)lli. < C'||5t;(t)||i2(||t;^(t)||i. + l) 



'V'^W 3 



+ Cexp{C\\\/v'^\\ ^ I )||po||i.o||(5t;||^i^2. 

Integrating this differential inequality we get 

\\Sv\\l^,. < C f\\6v{r)\\U\\v\T)\\l^ + l)dT 
Jo 

+ Cexp(C||Vt;l ^ | )||po||ioc / IMll^^^dT 

^t^p,i Jo 

< C f\\6v\\l^,,{\\v\T)\\l^ + l)dT 
Jo 

+ Cexp(C||Vt;l a )||po||L-t^ / \\Sv\\l^L'^dT. 

^l^P,i Jo 

It suffices now to use Gronwall inequality. 
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